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ABSTRACT 



We present the detection of deuterated molecular hydrogen (HD) in the remote Universe in a damped Lyman-a cloud at z abs = 2.418 
toward the quasar SDSS J143912. 04+1 1 1740.5. This is a unique system in which lU and CO molecules are also detected. The 
chemical enrichment of this gas derived from Znn and Sn is as high as in the Sun. We measure A , (HD)/2A , (H 2 ) = 1.5^S - * x 10 -5 , 
which is significantly higher than the same ratio measured in the Galaxy and close to the primordial D/H ratio estimated from the 
WMAP constraint on the baryonic matter density (f2 b ). This indicates a low astration factor of deuterium that contrasts with the 
unusually high chemical enrichment of the gas. This can be interpreted as the consequence of an intense infall of primordial gas 
onto the associated galaxy. Detection of HD molecules at high-z also opens the possibility to obtain an independent constraint on the 
cosmological-time variability of ji, the proton-to-electron mass ratio. 
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1. Introduction 

Deuterium is produced by primordial nucleosynthesis and is sub- 
sequently destroyed in stars. Therefore measurements of D/H 
from primordial gas provide important constraints on the bary- 
onic mat ter density Qb in the framework of Big-Bang cos- 
mology (Wag oner! 1 1973b . Measurements at different redshifts 
in turn provide important clues on t he sta r formation history 
dDaigne et all 120041: ISteigman etail 120071) . All the available 
D/H measurements at high-z are based on the determination of 
the iV(D°)/jV(H°) ratio in low-metallicity QSO absorption line 
systems. These measurements are difficult mainly because the 
velocity separation between D i and H i absorption lines is small 
(Avdi/hi ~ 80 kms -1 ) implying the lines are easily blended. An 
additional difficulty is the presence of the Lyman-o- forest, mak- 
ing hard to find the true continuum position and to discern be- 
tween D i absorption lines and intervening Ly-or forest lines. This 
explains why, despite more than a decade of efforts, only seven 
robust measurements of D/H at high-redshift hav e been per- 
formed dO'Meara et all 120061: iPettini et ail l2008al) . It happens 
that these measurements are consistent with the value [D/H] p 
= 2.55+0.10 x 10~ 5 derived from the Wilkinson Microwave 
Anisotropy Probe (WMAP) five year data dKomatsu et alll2008l) 
together with the r\ <-» [D/H] p conversion from Buries et all 



* Based on observations carried out at the European Southern 
Observatory, under programme 278.A-5062 with the Ultraviolet and 
Visual Echelle Spectrograph installed at the Very Large Telescope, unit 
Kueyen, on mount Paranal in Chile. 



(1200 11) . where 77 is the baryon-to-photon ratio. We consider this 
[D/H] p value as the primordial D abundance in this work. 

Damped Lyman-o' systems (DLAs) are absorbers with the 
highest neutral hydrogen column densities among Hi absorp- 
tion line systems (N(H°) <: 10 20 cm~ 2 ) and are thought 
to arise from the neutral inters tellar medium (ISM) of dis- 
tant galaxies (Wol fe et ail 120051) . Only a small fraction (-10- 
15%) of them show detectable amounts of H2 dLedoux et all 
120031: iNoterdaeme etail 120081) . Among the 14 high-z H 2 - 
bearing DLAs known to d ate, only two show H D absorption 
dVarshalovich et all l200lt ISrianand et all 120081) . The z abs = 
2.418 system toward SDSS J143912+1 1 1740 we present here is 
the only one where, in addition to H2 and HD, CO is detected. 
From the excitation of CO it is possible to derive in a straight- 
forward way an accurate estimate of the Cosmic Microwave 
Background Radiati on temperature at the redshift of the absorber 
(Srianand etal.,2008). Here we focus on the HD/H2 ratio and its 
implications for the star-formation history in the DLA galaxy. 



2. Observations and Analysis 

Observations were performed with the Ultraviolet and Visual 
Echelle Spectrograph of the Very Large Telescope at the 
European Southern Observatory on March 21-25, 2007 with to- 
tal exposure time on source exceeding 8 h. Both blue and red 
spectroscopic arms were used simultaneously using dichroic set- 
tings with central wavelengths of resp. 390 nm and 580 nm (or 
610 nm). The resulting wavelength coverage is 330-710 nm 



2 P. Noterdaeme et al.: HD at high-z: A low astration factor 




-2000 -1000 1000 2000 

Relative velocity (km s -1 ) 



Fig. 1. Voigt profile fits to the Lya (top panel) and Ly/3 (middle 
panel) lines. The bottom panel shows the absorption profile of 
O i to visualize the positions of individual H i components (ver- 
tical dotted lines). The solid red lines indicate the overall fit- 
model. Profiles are plotted on a velocity scale with origin set at 
Zabs - 2.41837. Blue profiles at v = kms -1 represent the main 
H° component and the associated uncertainty as a shaded region. 



with a small gap between 452 and 478 nm. The CCD pixels 
were binned 2x2 and the slit width adjusted to 1" matching 
the seeing conditions of ~ (X'9. This yields a resolving power 
of R = 50 000, as measured on the thorium-argon lines from 
the calibration lamp. The data were reduced using the UVES 
pipeline v 3.3.1 based on the ESO common pipeline library sys- 
tem. Accurate tracking of the object is achieved even in case 
of very low signal-to-noise ratio. Both the object and sky spec- 
tra are optimaly extracted and cosmic ray impacts and CCD de- 
fects are rejected iteratively. Wavelengths were rebinned to the 
vacuum-heliocentric rest frame and individual scientific expo- 
sures were co-added using a sliding window and weighting the 
signal by the total errors in each pixel. The dispersion around 
the wavelength calibration solution is 150 m s _1 . Standard Voigt- 
profile fitting methods were used for the analysis to determine 
column densities, redshifts, and Doppler parameters b. 

Hi absorption corresponding to the DLA is spread over 
~1000 km s . The velocity structure can be modeled using 
the asymmetry of the Ly-a line and the profile of the Ly-/3 
line (see Fig. [TJ. The derived structure is consistent with the 
velocity profile of the Oi absorption. The latter species is be- 
lieved to track H i closely t hrough charge exchange reactions for 
log N(H°) > 19.0 dVieeasl [19951) . The continuum is fitted at the 
same time as the absorption lines with a large-scale low-order 
spline. 

The column density in the main clump is strongly con- 
strained by the red damping wing and the gaps in the Ly-/3 
absorption. The best fit gives logN(H°) = 20.10!°^ for the 
main clump at z a b s — 2.41837. Note that a column density of 
log N(H°) = 20.2 would fill in the whole red profile (see Fig. 1). 
For a column density smaller than log N(H°) = 20, the shape of 
the profile can not be reproduced by any absorption located at 
-117kms _1 . 

The overall fit is represented in red on Fig.Q] while the con- 
tribution from the main clump only is shown by the blue profile 
and the associated uncertainties by the shaded region. Indicative 
column densities for other components are logA^(H°) = 17.9, 



deuterium in a solar metallicity DLA system at z = 2.418 
Table 1. Observed metal abundances 



Ion (X) 


logJV(X) 


[X/H]' 


S + 


15.27 ±0.06 


-0.03 ±0.12 


Zn + 


12.93 ± 0.04 


+0.16 ±0.11 


Si + 


14.80 ± 0.04 


-0.86 ±0.11 


Fe + 


14.28 ± 0.05 


-1.32 ±0.11 


N" 


> 15.71 


> -0.34 



1 with respect to solar abundances from Morton (2003). 



19.25, 19.20 and 19.20 at velocities v = -892, -594, -382, and 
-117 km s , relative to the main clump studied here. The com- 
plex velocity structure and/or insufficient signal-to-noise ratio 
in the blue prevent the determination of N(D°) in these indi- 
vidual components directly from D i lines. From fitting together 
various optically thin transitions of S n, Zn n, Si n and Fe n, we 
derive metallicities in the main clump (at v = kms -1 ; see 
Table [TJ. No ionisation correction is applied. The presence of 
strong neutral carbon and molecular lines in the main clum p, in- 
cluding easily photo-dissociated CO (Srian and et al.L l2008) indi- 
cates that the effect of ionisation on abundances should be neg- 
ligible. Indeed, in the abscence of neutral carbon, the ionisation 
correction for logjV(Hi ) = 20 .1 is smaller than 0.10 dex (see 
Fig. 23 of lPerouxet all 120071) . The S and Zn metallicities are 
solar and the relative depletion pattern (from Si and Fe) is typi- 
cal of what is seen in cold neutral ISM clouds in the Galaxy (see 
TableQ). 

3. Column densities of molecules and D/H ratio 

Absorption lines from more than one hundred H2 transitions 
from rotational levels J = up to J = 5 are detected in the 
main Hi component in six components spread over 50 km s _1 . 
Column densities in different J levels are obtained by simulta- 
neous fits and are especially well constrained by the presence 
of damping wings seen on the corresponding lines from the 
Lyman 2-0, 4-0, 5-0, 7-0, 8-0, 9-0, and 10-0 bands and from 
the Werner 0-0 and 1-0 bands. Examples are shown on Fig. [2] 
with^ 2 = 1.04, 0.99 and 0.98 for J = 0, 1 and 3, respectively. We 
measure a total column density of log N(H2) = 19.38 ± 0.10 and 
a molecular fraction, / = 2jV(H 2 )/[2jV(H 2 )+jV(H )] = 0.27+°J°. 
The uncertainty on NQI2) is dominated by the error from fitting 
the damping wings of low rotational level lines (J = and 1). 

Deuterated molecular hydrogen is detected in the first ro- 
tational level in three components associated to the strongest 
H 2 components. Five HD absorption lines (L3-0R0, L5-0R0, 
L7-0R0, L8-0R0 and W0-0R0) are clearly detected and were 
fitted simultaneously. The fits are shown on Fig. [2] Signal-to- 
noise ratios are -12 for L0-0R0, L3-0R0 and L5-0R0, and 
~6 for the remaining lines, x 1 values for the fit to these lines 
are 0.43 (L0-0R0), 1.04 (L3-0R0), 1.36 (L7-0R0), 1.07 (L8- 
0R0) and 1.49 (W0-0R0). Note that fittings were done inde- 
pendently by two of us (PN and RS) with two different tools 
(FITLYMAN and VPFIT), yielding same results. The strongest 
constraint comes from L3-0R0 transition which has good SNR 
and is optically thin (see Fig. 2). Note that L5-0R0 is blended 
with another absorption and is not included in the fit. Results of 
the fits are summarized in Table [2] Errors correspond the fit of 
both the lines and the continuum. The total HD column density 
is log MUD) = 14.87 ± 0.025. 

Although each HD component is associated to one of the 
H 2 components, the strong blending of the latter, especially in 
low rotational levels, does not allow for the determination of 
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Fig. 2. Molecular hydrogen absorption lines. The normalized 
flux is given on a velocity scale with the origin set at z a b s = 
2.41837. The Voigt profile fits as well as the residuals are also 
plotted. The blue dashed profile overplotted on HD L5-0 R0 cor- 
responds to the contribution from HD only. The grey regions for 
H2 L4-0 R0 and H2 L4-0 Rl represent the lcr uncertainty around 
the best fit. Short vertical lines mark the position of individual 
components. 

Table 2. Voigt profile fitting results for HD 
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Fig. 3. The HD/2H2 ratio vs the molecular fraction. The 
filled square is the new measurement at z a b s = 2.418 to- 
ward SDSS J143912+1 11740. The empty circles are fuse and 
Cope rnicus measurements in the Galactic ISM (lLacour et all 
I2005IK The grey region marks the D°/H° ratio in the Galactic 
disk (Linsky et al., 2006 ), the light gre y region marks this ratio in 
the Galactic halo dSavage et alll2007|) and t he vertically-dashed 
one t hat in the Local Bubble dMoos et all 120021: iLinskv et all 
2006). The oblique dashed lines stand for D°/H° measure- 
ments in_Jri^h-r^shift_quasar absorption-line systems (QSO 
ALS: lPettini et al.U2008al) . Finally, the solid line corresponds to 
the D/H ratio estima ted from the baryon/photon ratio (WMAP; 
iKomatsuet aiT.l2008l) . 





V 1 


b 


logA'(HD) 




(km s" 1 ) 


(kms- 1 ) 




2.41835 


-1.9 


2.9 ±0.8 


13.89 ± 0.08 


2.41851 


12.0 


5.0 ± 1.0 


14.57 ± 0.04 


2.41866 


25.3 


3.5 ± 0.3 


14.46 ± 0.03 



Velocity relative to z abs = 2.41837. 



7V(HD)/2/V(H2) in individual components. We therefore use the 
column densities integrated over the whole profile for both HD 
andH 2 and obtain V(HD)/2V(H 2 ) = 1.5!^xl0- 5 . InFig.[3] we 
compare the 7V(HD)/ 2A^(H2) ratio with s imilar measurements 
in the Galactic ISM ( lLacour et all 120051) . It is apparent that 
the V(HD)/27V(H2) ratio in the present system is an order of 
magnitude higher than the values for similar molecular fraction 
(and N(H2)) measured in the Galaxy. We also compare these re- 
sults to D°/H° measur ements in low metal licity clouds toward 
high-redshift quasars dPettini et all 1200 8a). the Galactic disk, 
the Galactic halo and the pri mordial D/H ratio e stimated from 
the five-year WMAP results (Ko matsu et alll2008h . 

From the molecular ratio, HD/2H2, we derive (D/H)dla > 
0.7 x 10~ 5 at the 95% confidence level. This corresponds to an as- 
tration factor / D = (D/H) p /(D/H) DLA < 3.6 where (D/H) p refers 
to the primordial abundance derived from WMAP. However, the 
true (D/H) ratio (resp. astration factor) is probably well above the 
derived lower limit (resp. well below the upper limit) for various 
reasons: 

First, the quoted HD/2H2 ratio could represent a lower limit 
on the actual ratio in individual components as H2 and HD 
are possibly not co-spatial. Indeed, the maximum HD column 
density does not ar ise in the strongest CO component (see 
ISrianand et alll2008h . suggesting a fraction of the molecular hy- 
drogen is not associated with HD. 

Second, deriving the deuterium abundance from the HD 
column density is difficult because of the complex chemistry 



(e.g. ICazaux et all 120081) and the sensitivity of the HD abun- 
dance on the particle density, cosmic ray density and UV field 
dLe Petit et all |2002|) . However FUSE and Copernicus observa- 
tions have shown that in the ISM of the Galaxy, the HD/H2 ra- 
tio increases with the molecular frac tion and that HD/2H 2 could 
trace D/H well only when / ~ 1 (lLacour et al ., 2005), when 
both HD and H2 are self-shielded from photo-dissociation. In 
diffuse gas, the HD optical depth is expected to be smaller than 
that of H2 -as the deuterium abun dance is low- and HD/2H2 
will p rovide a lower limit on D/H dLe Petit et all 120021; iLiszA 
2006 ). This is als o supported by recent FUSE observations by 
Isnowetal.l(l2008l) . 

Finally, the D/H ratio in the gas phase, as measured by 
HD/2H2, could itself be a lower l imit on the true abundance ratio 
if dep letion on dust is significant dProchaska et alll2005l:lDrainel 
2006). 



4. Discussion 

Evidence for infall? 

The high D/H ratio inferred above indicates that astration of 
deuterium is low eventhough metallicity is solar. This situation 



primordial material ( 


Steieman et al!l2007llRomano et alll2006; 


Prodanovic & Fields 


, 2008b. 



If we use the sulfur abundance as a proxy for that of oxygen, 
we note that the [N/O] ratio in the main H 1 component is con- 
sistent with the ratio expected for secondary nitrogen production 
at solar metallicity dCenturion et all 120031: IPetitiean et alll2008l: 
iPettini et all l2008bl) . In addition, the [S/Zn] ratio is consistent 
with the solar value and does not indicate any a-enhancement. 

The inferred lower limit on D/H, the near solar value of 
[a/Fe] and [N/a] rule out rapid star formation that is generally 
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invoked to explain high chemical enrichment in elliptical galax- 
ies. 

Speculating slightly further, we note that the metal profile is 
spread over at least 800 kms -1 (see Fig.[TJ. This, together with 
the solar metallicity, is co nsistent with the g as being associated 
with a deep potential well dLedoux et all l2006h . 

The properties of the gas in the present system are similar to 
those of the ISM in our Galaxy. But these properties are reached 
on a time-scale five times smaller for the DLA galaxy than for 
the Milky Way, as the age of the Universe at z = 2.42 is 20% of 
its pr esent age adopting the most recent cosmological parame- 
ters (iKomatsu et alll2008l) . This must imply that the system has 
undergone continuous star-formation and infall. 

If one assumes negligible production of deuterium by cosmic 
rays and that this element is completely destroyed in the mate- 
rial that goes through stars then the gas infall rate M[ n should 
on average be of the order of the star formation rate Msfr in 
order to replenish deuterium. Recently several observational ev- 
idences have been published for cold gas accretion onto massive 
galaxies at high-z (|Weid inger et all l2005t iNilsson et all l2006t 
IDiikstra et ail 120061) and low-z dFraternali & Binnevi l2008O n 
addition, numerical simulations suggest that at z — 2 - 3 the ac- 
cretion of cold material fr om the IGM domin ates for halos with 
masses <~ 3 x 10 11 M Q (iKeres et all 120051) . Interestingly, the 
inference that M; n ~ Msfr is also required t o understan d the 
properties of z = 2 - 3 Lyman-break galaxies dErbl 120081) . Our 
observations reinforce this important finding. 

Constraint on the variability of /j 

The detection of several HD transitions should make it possible 
to test the time variation of // = m p /m e , the proton-to-electron 
mass ratio. H2 transitions at high redshift have been used (e .g., 
Ilvanchik et afll2005hlRernhold et alll2006tlKing et all 12008b to 
probe the variability of fi. This is done by measuring the relative 
position of the lines around the overall redshift of the absorbing 
cloud, 

£ = fo - Zabs)/(1 + Zabs) = — Kj , 

where z, = - 1 is the observed redshift of line i and 

Kj — d In A® /din fi is the sensitivity on fi calculated for each tran- 
sition. z a bs is taken for each component as the weighted mean 
redshift from the different transitions. As these measurements 
may involve various unknown systematics, it is important to use 
different sets of lines and different techniques. Sensitivity coef- 
ficients and accurate wa velengths for HP tr ansitions have been 
published very recently (llvanov et alll2008l) . However, the SNR 
of our data at the position of the HD absorption lines prevents 
us to perform a measurement at the level of current studies. 
Additional data on this quasar are needed to be able to derive 
an independent constraint on the variability of fi. 

5. Conclusion 

We reported the detection of HD molecules at z a bs = 2.418 to- 
ward SDSSJ143912+1 11740, following a careful selection of 
quasars in the SDSS database and intensive observations with 
UVES at the Very Large Telescope. The system presents charac- 
teristics very similar to what is observed in the solar neighbour- 
hood. We find HD/2H2 = 1.5 x 10~ 5 , consistent with an astration 
factor of deuterium less than 1 .7 which is contrasting with the 
high chemical enrichment. This is best explained by a scenario in 



which the gas that goes through star formation is replenished by 
continuous infall of ambient primordial gas. Similar results arise 
from re cent numerical simulations a nd semi-analytical mod- 
els (e.g. lKeres et all 120051; lErbll2008l). Interestingly, d ynamical 
studies of nearby galaxies {fraternali & Binnevi 120081) as well 
as int e rpretation of Ly-q blobs at high redshift Jyy eidinger et all 
120051: INilsson et all l200d IDiikstra et all 120061) provide inde- 
pendent observational evidences for accretion of gas onto mas- 
sive galaxies. 

We finally stress the importance of detecting similar systems 
to probe the time-variation of the proton-to-electron mass ra- 
tio from HD lines. Although such a independent test would be 
welcome to characterize possible unknown systematics, signif- 
icantly higher signal-to-noise ratio is required to obtain limits 
comparable to those obtained with other techniques. 
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